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INTRODUCTION 


when Einstein first introduced the notion of the geodesic of 
@ particle in the General Theory of Relativity (1), he did so 
through the principle of equivalence, postulating that the motion 
of a particle under the influence of a gravitational field must be 
indistinguishable locally from inertial motion, The geodesic 
equations contain the Christoffel symbols of the second kind, field 
variables which effectively take the place of the force of classical 
theory. The question arises, whether these equations can properly 
be called equations of motion, if the Christoffel symbols are sin- 
gular at the location of the particles for this is the case with the 
only solution of the field equations in empty space which is static, 
has spherical symmetry, and which goes over into the flat metric at 
infinity, ise. the Schwartzschild solution, 


In 1938, Binstein, Infeld, and lioffman(2,3) collaborated in a 
series of papers in which the necessity of considering the above 
question was obviated, Thay succeeded in deriving the equations 
of motion of a particle in a gravitational field directly from the 
field equations, without bringing in the idea of a geodtsic at all. 


In this paper, an investigation is made of a more general set of 
field equations, derived from a Hamiltonian principle using an integ= 
rand slightly more complicated than that which gives the accepted 


mS 


equations. This set, on the basis of the principle of simplicity, 
has not been exhaustively examined previously, for it involves fourth 
derivatives of the metric tensor. However, it will be shown that 
these new equations reduce to the original ones as a certain para- 
meter tends to sere, and moreover give, for the static case, a solu- 
tion which is regular at the location of the particle, So we find 
that the equations of motion, which are geodesics, follow directly 
from these field equations. 
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THEORY 


We are looking for a aymetric second-order tensor, pel 
depending only on the metric tensor and its derivatives, which we 
can identify with the stress-energy tensor 7°” 4 The principle 
of conservation of sxorgrcamentan stétes, Nowerer, tes 77", J= 0 3 
So our second~order tensor must satiafy 1.4 = 0 whore 5 sige 
nifies covariant differentiation. ( the notation , is used to denote 
regular differentiation ). 


Following Sehridinger (4), wo show in Appendix A that if 77 
is any scalar density that depends only on the Aik and their der« 
ivatives with respect to the coordinates up to any finite order, 
then the invariant divergence vanishes identically for the tensor 
density that is constituted by the Hamiltonian derivatives of )7/ 4 
feo, ot T= [74 dx* +» Compute 3 C and perforn 
the partial integrations needed to put it inte the form 
JL=S(M/09,4) 00,4 det « tat 176° = S10 / 14. » Then 
mR = () por, (11/5); = O « So wo take yan" equal 
to mF and identify this with 7". 


We shall wish to consider, as the integrand in our Haniltenian 
principle, the product of /-7 with F , where is a scalar 
depending only on gi* and X:k as independent variables (f;/, 
is the Binstein-tensor ). Using the notation, 2 ‘* = OFT Ritk 9 
we find, (Appendix B) | 


“lest 


Ln = t (4 p' ‘+ gh a gp a aig} sh fe 

+g" bond cps (1) 
therefore, following our previous analysis, /.°* = se agent where 
‘ih oe ig the expression in parenthesis, and ‘/ is a constant to 
be deternined, 


We are now ready to consider a definite form for / « Instead 
of taking F = RV-9 4, a8 Einstein did, we consider the next 
simplest expression, - = R+aR® +bRy RH » whore 
R is the curvature scalar, and « , 6 are arbitrary parancters, 
Remembering that R= 9S Ry where Rand G4 
are our independent variables, and that Keg ee » 
we have 
IR = 9% JR = akg" TRyRY _ 5 aih 
DRik PRik ate 


k k é 
ae’ gh SR __apr® SRyRI__ gh aH 
J Paik Ty “aq 


Substituting into equation (1) we get (Appendix C), ; 
2 b Maes A jo 2 ——t 
Bog (R)+ & Ape (R*) +b Ape (Ry RI)= 195 Jno 


Pg (2) 


= 


where Bpg (R) =-2 (Rpg ~ > A pg R) 
Deg (a2) = 2 (g/g “64 getght_g 99 gg") Rae 


is argis — yYRRER 


Ag (RyuRM)= aes +f yeg"* 
- 4, G9 aoe ML "A )Rydge 


Rea Sei ali heintie 


Now assume that the gravitational fields (is@., the deviations 
of the actual metric from the flat metric) encountered in nature 
are So weak that the non-linear character of the field equations 
leads only to secondary effects, So we apply the linear approxi- 
mations 

jie 
where A is the parameter of expansion and a small constant, Epo 


= Eon tAheot A hpm ts 


is the flat mgrie /-/ © © o where we have taken ¢*- / « 
o-/ Oo 
D o-! 0 
o O O+/ 

To first order in ) ina aystenwhere J- = O 4, we 


Calculate (Appendix D), 


i rg R= -AE “SV 4g, af 
Qpg (R* =~ 216 9g é* “2 SVarele t2rETY ar pg 
a ta Vice AE a a 
Og (Ry RY) = Mle“ Tagapte— pz Le Yeap oe 


x E°SY arpa] 


Parag ' 
and “7 ) 41s defined by hinge Vur-¥ L..) 
ey 
T= k= eh) 


Substituting inte equation (2), we get 
~Aé€ <7) af + BnA (E*F 2, appa Ege EF V apse) 


tbhrle LY appgy ~€ ge ti “LY apele 
ewe eo, 2 pce) = 1" ps pg (3) 


If we take a ates aie 
lead to equations which can be solved conveniently, and to which 
regular solutions can be chogen, So, purely for reasone of con- 
venience, we choose b=24 4 


tes 


This leads to the equations 

o picé r a 
-2ES V9, ap 2a NEWEST Dog, ape = 4 Tog - (1) 
Equations () constitute our field equations. 


Let us now verify that one can choose a solution of these equatiom 
in empty space which is static, has spherical symmetry, which goes 
over into the flat metric at infinity, and moreover is regular at the 
location of the particle, 

Inemty space 7,,= 0 yard P= Fyy 0 the only 
component of importance, Putting 22-=-//h*, wire k>o , 
the field equations fob the static case become 
viv G-k 7 g= O 


cA aprepe > coordinates, J 
Fle D-igM)- 0 
The general solution to this equation is 
d-# ety BoM) S43; 


For () > aH as 2 —>.4 » we take B= D=0, and then 
eo 
b-=fe™“+ 
wo ye be ae Oe 
= le “RAL 
‘inins 0s eine li? O 8 
w@ Q> os k-—>  » our equations reduce to Binstein's 
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equations, and approaches -C /1 y the Schwartaschild solution, 
ifwelet C= 24 / where 17 4s the mass of the point particle 
We are considering, and © is a wiiversal constant which has the 
wee G=6.4x10° tues i 

Since we can always find a coordinate system in which the pare 
thele is at rest, will be regular at the origin in the dynamic 
case as well, It remains to note that the derivatives of p are 
not analytic at the origin, for the;depend on the direction of app= 
Toath. It is possible to circumvent this problem by averaging over 
all directions, but in the final analysis, we would need to consider 
field equations containing fourth and sixth derivatives of the met- 
ric tensor, so as to be able to apply boundary conditions to the 
first derivative of p . 


It is possible to show that the geodesic equations of a pare 
ticle follow directly from the fact that the invariant divergence 
of the stress-energy tensor is BEM Lees 9 Bigs iso . 
The procedure is as follows, 
We note that (Appendix 8) 


4 at | 
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Wow (Appendix F) a 4 A ee x 

-gT? = on fo EC Ox 2ads 
Therefore 


D (FaT VT) = mf dz d2?d Si _xe@)ds 
2, FG ) mf SE EE, v2) 


aQus 
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fae de? a SA(x- 2O\ ds 
=o). Ge Bee 


‘ -m f a3 of 54 (x2 (s)) os 
pany wapicsapeehcete agra ene 
m [ox "Ede ola = m[ Teed e* oe wee 

: m{o (x-2)fe- rds eas dz! 

=m) (x-2)ds dX 

” ts ot a.*™ 
evaluated at 2 '=7 , 
A 4 
tare FF AhTP = is mf 2“ de lz” 54 exo) ols 
- 3 Ymf dada plea lea oft 
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evaluated at ads 


rae ae 

. z 
wox«*xZ ee ae xe ae » then, 
— is not zero, and the deltasfunction is infinite, 


rae + 15 & hae ees O 
) a (2 


weit e, ae 
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if the Christoffel sybols are regular at the location of the 
particle, we can say immediately that these equations are the equae 
tions of motion of the particle. We cannot take this step in the 
original theory, however, for this in not true in general, But with 
our new field equations, we know that the Christoffel. symbols are 
regular, Our object now is to solve these field equations to first 
order in “1 , with the solution in integral form, and to see what 
happens a3 a—> 0» 

We want the solutions to the equations 
- AES 9 pg, af —~20A E** °F Nog apne eT ea 
Taking 2/°= 1/20 and -" /2a =A 5 this becomes 
(a*+0/*Q) y? . AT 
were [J] = 1%42-V~ ; 
We are permitted to raise the indices uith the yo as the error 
in go doing is of the secomd order in \ , The Green's function 
C, (%,%¢') Sor thie equation is defined by 

a FT et sie ere e} 


ie olp, dps dpe ~7* 
G aie oe ary ell p< prP - Qr p-p 


as can be seen by direct substitution, This can be factored to give 


OS) = Bay [A is scale ee 


will 


al p eth (XA) 

- {ff ftage tema 

nce nbbaitieltsdililé-d tanttee ealuieths ech cell uneaiiibes 
CUES Senate 


mere x~) 
ab ff ere pu ~ 


aa 


seeds d, (Xu) Uy > Ua 
= Ke lUgl< UA 
O 


Un<-Un 

, = / 3 j Yo 
where ee a, ? Un= = (- uu)? 
ma (>= Uy —-Urn « 


Substituting, wo get for our retarded Green's function Gr (Xp xe) 


= tf (25) Beal, Yar) 


Ay PUA 
O ~— O [Ug <td 
O O Uy <A 
it Gal \ tae 
O lua | <Ua 


Wi) =A [45 G(x, IT () 


and 
HT bg )= far 2 O10G- aoe 
Therefore 
0G) = Am fob ot a lla 8 39S l- 22) 


= Aw iC 7 2A () 2" (s) 
where in the last expression (y= X¢~ 2 « 
We now, for convenience, introduce another constant, i » defined 
ty bra 2cSsia, mem Az~7Rad=-7/b*A 


gy tel ay HE wet 


(=a 
This is the retarded solution.in integral form that wo wore 

looking fors We know that as a-> 0, equation (4) approaches that 

of Binstein if we take // = (bm G (Appendix G)_ Let us now 

see if we arrive at his solution if we let 2 approach zero in our 

solution, In this case, /-> © as welll, and in onder to havo a 

finite answer, we must let m approach zero with b such that m/b is 

a constant, m', Aloo, as b-70, -1, (uls)/b) /u le) becomes a 

Dirac delta~function and the integral can be solved exactly to give 

1“) = tw ete ye ie) 

A K isa ls (Appendix H). 
where {| is the proper time of the retarded point, ise. the point 
on the world line euch that (/ i= © and lly > CO « In the 
Cutan tie shah, alee On Gemmmaihin Wh th & Seiatien ot S is 
not under an integral with respect to SS , then it is taken at 
f= So « K 4s defined by 


=-13- 


k= ~ Ua du*/da. = (Xe-2y)H* = Ux 
where a dot above a function of 5S denotes differentiation with rese 
pect to S pend 17 = Z 76 
fo first order in A , V-9 = 1 ARE A pe » and so 
to this order, Y"0G) = -(1/A Mn’ 7R) uAy (7) 
fo fix my Look at the ly) component of 004) 4m the static 
case, and compare with the expression given in the original theory, 
eee 


; a \ eh we av [mia we “2)W_ —dGm F 
AV yyl) 46) he eth - 4G lA-A1l [A- 2 
uf 
Sey: (0,0,0,1) »mma K=Uy=U 6 see, 
mince UU“ =O» Uy = SZ U; gt y2p3)e We aloo 
note that liy~ X,,~ 2, (S) 4s powitive, aince we are considering 
values at. the proper time of the retarded point, Then we have 


taph nw EEE / ~ > 

» gs om i © is aga er - — iG PEt. 
A} U4 vt Xe ~ 

Comparing these two expressions, it ii imas ‘= ms (1T 


Wo will now exhibit the Christoffel symbols derived from 

is ; 
Why Bly Ne Letting 6 = 2am » we geti/appendixc 1) 

= fo ale U4 OT U, V5 vFe Uy. vue vf 

v5 ula, yuk Vy Ye VET 
I-k’ [ 4-4-f f pu: vu 

eh ee i 

ae re) [é.,u* Js =, Na] 


+ Aa 
~t£ 55 [ Gus 204 y. 7} " 


where X “= 4.07 
As a check on this formula, let us again consider the case of 
a stationary particle, Equation (8) reduces to 


ae 


A lrligu, tr rly yee!) 


+ ++(é, Us— Uy o.. fu, d,") 
Due 


- FEav eT yd!)-vy x) 


Tals gives, where 7 > 4.4) % 4 » 


LP bal til bel iif Se 
GP LEE Meh lhe Lo PLO 


Also, under the above conditions, it is easy to show that 


Y. = — 228 ; 
V4 rd, f A (U4)? 
g -2 HU = AL Y4d, + 


Uj” OF! 
In our theory, served as an expansion parameter, Neverthe- 


less, we may now choose a definite value for it, and we take it equal 
to unity. Then the Christoffel symbols check with the values cal- 
culated by Minstein, Infeld, and Hoffman (see, Gegey Bergmann p.230). 


It is now possible to consider the case of two interacting point 
particles. Let =, = label the coordinates of particles 1,2 res= 
pectively, and let % label. the coordinates of the field point. 


-15= 


Thea, if , 5 2. de the world-Line parameters of the two particles, 
we have, writing dow the geodesic equations of each of then, 


a | 
d 2 a { fixes) 12,“ da? _ Oo 
Aa > ‘aia bu cae clo. “be, 
Pare 


des +4 hy eta Az — o 
ng Lf 0" Chay “olor 

Now let 2% > =. in the first equation, and let x > =, dn 
the second, These then represent the equations of motion of the two 
particles. 
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SUMMARY 


it was shown that the invariant divergesce venishes identically 
for the tensor density that is comptituted by the Hamiltonian deriv~ 
etives of any scalar dentity that depends only on the metrie tensor 
and its dovivativas with reapect to the coordinates uw to any finite 
order, This tensor density could then be identified with the product 
of the atress-omergy tansor with the square root of the detenainant 
of the metric tensor, this identification to be considered as a 
possible set of field equations was derived which reduces to Einstein’ 
equations as a cortein parameter “a" tends t sero. 

The linear approximation was applied, ani a solution was found 
for the static case which wes regular at the location of the particle. 
This moant that the Christoffel aymbols were also regular there, It 
vag then shown that the geodesie equations of a particle follow 
directly frou the vanishing of the invariant divergence of the stress- 
energy tengor’y oF, We can say, dlrectly from the field equations 
Sheuselves, These equations contain the Christoffel syubols, which 
wore showm to be regular at the location of the particle; so the 
geodesic equations can now be called equations of motion. As the 
parameter "a" tends to zero, the Christoffel symbols approach sin- 
gularity in a regular way; ‘this could be called a regulator technique. 


In the latter part of the papes, the time-dependent field 
equations were solved in the linear approximation with the solution 


oe 


in integral form, Then, letting "a" aypreach sere, the integration 
was Carried out, and the solution was shown to approach that of 
Einstein for the static case, Also, the Christoffel symbols were 
exhibited amlicitly, end the case of a stationary particle was 
given special consideration, 


Ia con¢lusion, the new field equations lead te a Mewmian 
Geometry satiatying the principle of equivalence in which the equa- 
tions of motion follow from the field equations and are geodesics. 


APPENDIX A 


To show that if y 

ss {oa dex 

where 7 is any scalar density that depends only on the Fik 

and thelr derivatives up to any finite order, and we compute J 
and perform the partial integrations needed to put £t into the form 


= pe 4 
‘sh f d Fik 19 ik dx P (Ae1) 


then ‘ete O where m*. Tmt /d9 ° 


since [L te the difference between tuo invariente, it itecif 
is an invariant. Then me? 9:4 is a scalar density, But 19. » 
being the difference of two symmetric second rank tensors, is itself 
a symmetric second rank tensor. te 71" gy uw tented tinent 
rank, contravariant tensor density. 

Now we bring about a variation in 4; by changing the frane, 
toon wore X= %0 + AM (0) + A> Ue (LE) 
where A 4s an expansion paraneter, This expansion approaches iden« 
tity as /|—> 0. The integral then becomes 


L= JH Gen 14), 0) ol 

7 4 ef aX (X-) 
where 7’/ 419 the same function of the 9 bu (XL )= ox OF, / . 
and their derivatives wits respect to the > as it was before of 


~L9= 


the Fim») and thet derivatives with respect to the X » The 
limite of integration are the same if we ask that the transformation 
shall approach to identity at the boundary for any 1 . The only 
formal change is that the argunont is now not 92h) wa 95), 0) 
Expanding with respect to A , 

Wik 1-9; ke) = AL fa rad + fim He, ts oak é7 + 0A*), 
all functions written with the »° « Since this must vanish at the 
boundary for any A 4 use equation (4.1) with first order terns 

in A for ¢. :h » Dvopping superfluous dashes we get 


IZ = [770° (4 2H + 9m dtm + 244d )dr"= 0. 


Ox~XkR OXy 
Integrating by parte, ' 
De J A= kde wh ta ae 
OX; OD Whe Sa 


j k th 4 
= f(-2 Pe + A Bee \ dant. 


Now j, is arbitrary, hence 
2, oe ae 29:k oo) 
OX. yen 

Using the symetry of m*, 


ik. ¢ n ck oe ’ ™ 
er ee 


= m: Shim f. 
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3 kind, Thus 
whore fait is the Christoffel symbol of the second 
' Be re 
ute m*, ifs a oO, 
OX 


fou ete 
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APPENDIX B 


To show that if we take 77/ = (- ( ik Rek4 in our 
Hamiltonian principle, we arrive at the identity 

r hi xe, ke Xk LE ‘ = 
lgig’ OG GE hs age LE? 
where Fa = OF 

: dRik 


Consider [> = Sr vz dx”. We comute oT « 
OS IG + (USE pl) 9 
+ fr Plig lant 


We shall call these the first, second, and third integrals respec- 
tively. Use is made of the following relations: 
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following the rule for covariant differentiation of a tensor density 
(Schrodinger (1), Pe33)e Using the symetry in < and ( in [A}s 
and oti foc and ME ears cima 
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these results, 


d= dL = [f-€ L£(9pER!, gytek, ve ) 
bers eas igtdf bac - 


aes oS: ne * « Se, 
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APPENDIX 6 
To show that taking F = R pak + ba RY baad tae 
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Substituting into equation (1), we have 
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APPENDIX D 


To show that in the linear approximation, in a system in which 
Tn = Gp 
wit 
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We can find (see Bergmann (6), p.183) a coordinate system in 
wieh = ETD. = 0 
8 Radste = DEST 3 apueLéy Yast] SEM ante 
in a system in which v= Oe 


a Ay 
We have insidentally show that to this onder, eKyR » 
a ae » and Ry Ry all vanish, 
So, in the linear approximation, 
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To show that in the static case, the hh component of 7 
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To show that the Christoffel aymbols derived from expression (7) 
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